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We  present  a  quartz  erystal  mierobalanee  (QCM)  based  system  for  high 
sensitivity  differential  sputter  yield  measurements  of  different  target  materials 
due  to  ion  bombardment.  The  differential  sputter  yields  ean  be  Integrated  to 
find  total  yields.  Possible  Ion  beam  eonditions  Inelude  Ion  energies  In  the  range 
of  30-350  eV  and  ineldenee  angles  of  0-70“  from  normal.  A  four-grid  Ion  opties 
system  Is  used  to  aehleve  a  eollimated  Ion  beam  at  low  energy  (<100  eV)  and  a 
two  grid  Ion  opties  is  used  for  higher  energies  (up  to  750  eV).  A  eomplementary 
weight  loss  approaeh  Is  also  used  to  measure  total  sputter  yields.  Validation 
experiments  are  presented  that  eonfirm  high  sensitivity  and  aeeuraey  of  sputter 
yield  measurements. 


I.  Introduction 

Sputtering,  the  erosion  of  a  target  material  under  particle  bombardment,  first  discovered  in  the 
middle  of  the  nineteenth  century,  is  still  an  area  of  active  research’.  Much  of  the  experimental  work 
aims  to  measure  the  total  sputter  yield;  however,  the  angular  distribution  of  sputtered  particles,  or 
differential  sputter  yield,  is  also  of  great  interest  from  both  fundamental  and  applied  points  of  view. 
Differential  sputter  yield  measurements  are  important  for  better  theoretical  understanding  of  atomic- 
scale  collision  and  emission  processes’  and  are  also  of  practical  interest  in  areas  such  as  sputter 
deposition^  and  space  electric  propulsion  (EP)^'’"’.  In  the  former,  knowledge  of  the  angular  distribution 
of  sputtered  material  is  essential  for  optimization  of  manufacturing  processes.  In  the  latter,  angular 
distributions  are  necessary  to  analyze  re -deposition  of  potentially  damaging  material  sputtered  by 
energetic  ions  from  plasma  thrusters,  for  example  deposition  of  sputter  products  onto  solar  panels  or 
thermal  control  surfaces  can  compromise  spacecraft  operation.  Furthermore,  the  sputtering  in  many  EP 
thrusters  is  predominantly  due  to  low  energy  ions  (-30-100  eV)’^’’"’  for  which  the  corresponding  sputter 
yields  are  relatively  low  (<~1  atom/ion  in  atom  based  units).  Therefore,  the  ability  to  make  sensitive 
differential  sputter  yield  measurements  at  low  energies  is  critical. 

Total  sputter  yields  can  be  measured  using  weight  loss’^,  profilometry’^,  quartz  crystal 
mierobalanee  (QCM)'’'^’”'’^,  and  optical  emission  spectroscopy  (OES)^”.  The  high  sensitivity  of  the 
OES  and  QCM  techniques  allow  sputter  yield  measurements  at  low  ion  energies  near  the  sputtering 
threshold.  Total  sputter  yield  measurements  of  argon  ions  on  gold  were  measured  using  the  QCM 
technique  for  ion  energies  of  10-100  eV’”  and  total  sputter  yields  of  xenon  ions  on  molybdenum, 
titanium,  beryllium  and  carbon  were  measured  at  energies  as  low  as  10  eV  using  OES^”. 

Measurement  of  differential  sputter  yields,  especially  for  low  ion  energies,  is  more  challenging. 
Differential  sputter  yields  are  most  commonly  measured  using  collector  surfaces’^'^’'^’’  that  are  used  to 
collect  the  sputter  deposited  material,  which  is  later  analyzed  with  different  analytical  techniques,  such 
as  photometry^^,  profilometry^^,  Rutherford  backscattering^^'^*,  secondary  ion  mass  spectrometry^”’^^, 
electron  probe  microanalysis^”,  radioactive  tracer  measurements^^  etc.  The  use  of  collector  plates 
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precludes  in  situ  measurements.  Other  methods,  such  as  quartz  crystal  microbalance^"*’’°"’^’^’"^^’,  mass 
spectrometry^"^'^^,  laser -induced  fluorescence  (LIF)”,  and  cavity  ring-down  spectroscopy  (CRDS) 
also  used  for  differential  sputter  yield  measurements,  are  free  from  this  disadvantage.  A  disadvantage 
of  the  mass  spectrometry  technique  is  its  significant  geometric  limitations.  Mass  spectrometers  are 
usually  mounted  at  a  fixed  location  relative  to  the  vacuum  chamber,  which  limits  the  ability  to  measure 
full  angular  distribution  of  the  sputtered  material"^^.  The  optical  techniques  of  LIF  and  CRDS  are  both 
attractive  owing  to  their  non-intrusive  nature,  but  have  limitations,  for  example  the  need  to  match  the 
wavelength  of  the  laser  source  to  the  species  being  studies.  QCM  based  studies  are  further  discussed 
below.  Unlike  total  sputter  yields  which  have  been  measured  at  ion  energies  close  to  the  sputter 
threshold  for  some  ion  -  target  combinations,  there  is  a  paucity  of  differential  sputter  yield  data  for 
low  energy  ions  (<~100  eV).  The  majority  of  studies  have  been  performed  in  the  ion  energy  range  of 
about  one  thousand  electron  volts  and  above  though  some  studies  have  also  been  performed  at 

the  lower  energies  i.e.  50-1000  eV 

In  general,  experimental  results  of  differential  sputter  yields  obtained  at  energies  above  ~1  keV  are 
in  good  agreement  with  prediction  of  Sigmund’s  theory,  i.e.  sputter  yields  with  no  azimuthal 
dependence  and  simple  cosja)'’  dependence  where  a  is  the  polar  ejection  angle  and  w  is  a  constant 
greater  than  unity  for  over  cosine-profiles  and  less  than  unity  for  under-cosine  profiles.  Approximating 
such  profiles  with  a  simple  diffuse  form  proportional  to  cosja)  may  be  adequate  in  some  cases.  On  the 
other  hand,  for  energies  below  ~300  eV  systematic  deviations  from  Sigmund’s  theory  have  been 
observed*’"^^.  These  profiles  (for  non-normal  incidence)  tend  to  exhibit  complex  shapes  that  additionally 
depend  on  the  azimuthal  angle,  typically  with  increased  sputtering  in  the  forward  direction  and 
reduced  sputtering  in  the  backward  direction.  (Full  geometric  definitions  can  be  found  in  the 
aforementioned  references  but,  approximately  speaking,  the  backward  direction  is  back  towards  the 
incident  ion  direction  while  the  forward  direction  is  that  corresponding  to  specular  reflection  of  the 
ion.)  In  such  cases  the  use  of  a  diffuse  profile  does  not  provide  adequate  description  and  full 
measurement  is  critical. 

As  mentioned  above,  QCM  based  systems  can  be  used  for  both  total  and  differential  sputter  yield 
measurements.  When  the  QCM  is  used  for  total  sputter  yield  measurements,  its  surface  is  coated  with 
the  target  material,  and  then  bombarded  by  the  ions  to  measure  the  mass  removal  rate,  which  is 
proportional  to  the  sputter  yield.  In  contrast,  for  differential  sputter  yield  measurements,  the  QCM  is 
used  in  deposition  mode  meaning  that  it  acts  as  a  collector  surface  and  is  coated  by  particles  sputtered 
from  the  target.  Instead  of  a  single  measurement  of  mass  removal  rate,  the  mass  deposition  rate  is 
measured  at  multiple  sputtering  (ejection)  angles  by  scanning  the  QCM  above  the  sputtering  target. 
Only  a  limited  solid  angle  is  collected  at  each  location  (typically  about  1/1 00th  of  hemisphere)  so  that 
higher  sensitivity  in  the  rate  of  mass  change  is  required  for  given  ion/target  conditions.  The  mass 
deposition  rates  allow  calculation  of  the  differential  sputter  yield  at  the  corresponding  angle,  and  a 
series  of  such  measurements  can  be  used  to  describe  the  full  differential  sputter  yield  profile.  The 
differential  sputter  yield  can,  of  course,  be  integrated  to  give  the  total  sputter  yield.  Compared  to  other 
experimental  methods  used  to  measure  differential  sputter  yield  measurements,  the  QCM  approach  has 
several  advantages  including  the  possibility  to  perform  in  situ  measurements,  high  cost  efficiency,  and 
relative  simplicity.  In  our  apparatus,  in  order  to  measure  differential  sputter  yields  with  reasonable 
angular  resolution  at  low  ion  energies  as  is  needed  for  EP,  we  have  focused  on  high  sensitivity.  The 
system  can  measure  mass  deposition  rates  as  low  as  ~10'^  g/C  allowing,  for  example,  differential 
sputter  yields  for  ion  energies  as  low  as  60  eV  for  Xe  ions  on  BN  which,  as  described  above,  is  critical 
to  EP  applications.  In  contrast,  past  work  using  QCMs  for  differential  sputter  yield  measurements  has 
been  primarily  at  higher  energies*’^^’^^. 

Sputter  yield  measurements  at  low  ion  energies  require  appropriate  ion  sources.  For  example,  past 
work  using  QCM  sensors  employed  a  magnetron  as  an  ion  source^^’^^.  The  advantage  of  the  ion  source 
described  in  this  contribution  compared  to  the  magnetrons  from  the  studies  mentioned  above  is  that  the 
ion  source  allows  one  to  obtain  nearly  monoenergetic  ion  beams  over  a  wide  energy  range  (from  near¬ 
threshold  energy  to  hundreds  of  electron  volts)  and  with  different  ion  incidence  angles  relative  to  the 
target  (while  magnetron  sputtering  is  always  at  normal  incidence).  The  ion  beam  also  has  very  well- 
defined  properties,  including  low  beam  divergence,  narrow  energy  spread,  and  low  doubly  charged  ion 
fraction,  which  are  essential  for  high  accuracy  differential  sputter  yield  measurements. 

In  this  contribution  we  describe  a  quartz  crystal  microbalance  (QCM)  based  experimental  system 
for  high  sensitivity  differential  sputter  yield  measurements  in  the  ion  energy  range  of  30  to  750  eV.  The 
measurement  approach  uses  a  combination  of  weight-loss  and  QCM  deposition  sensor  techniques'’'*’’®' 

In  Section  II  we  discuss  the  experimental  apparatus,  while  Section  III  presents  measurement  and 
data  analysis  procedures.  Validation  and  demonstrative  measurements  illustrating  the  sensitivity  and 
the  accuracy  of  the  experimental  system  are  presented  in  Section  IV.  These  results  include  total  and 
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differential  sputter  yields  of  gold  under  bombardment  by  xenon  ions  at  ion  energies  range  of  30-350  eV 
and  at  incidence  angles  0-67.5°.  Conclusions  are  given  in  Section  V. 


11.  Experimental  apparatus 


A.  Description  of  Sputter  Measurement  System 

A  quartz  crystal  microbalance  (QCM)  deposition  monitor  is  used  to  measure  differential  sputter 
yields.  Total  sputter  yields  can  be  found  by  integrating  the  differential  sputter  yield  profde  and  are  also 
directly  measured  using  a  weight  loss  approach.  A  schematic  of  the  experimental  apparatus  is  shown  in 
Fig.  1.  The  main  elements  of  the  system  have  been  previously  descrihed*'*’  In  this  subsection  we 
give  an  overview  of  its  essential  features,  while  the  following  subsections  detail  specific  aspects  and 
modifications.  Improvements  relative  to  our  past  work*  that  enable  measurements  at  low  ion  energies 
and  yields  include  use  of  a  new  four-grid  ion  source,  higher  accuracy  temperature  control,  and  use  of  a 
new  QCM  crystal  and  controller.  The  new  source  allows  a  collimated  monoenergetic  beam  at  the 
needed  ion  energies  while  the  other  modifications  reduce  measurement  noise  thereby  allowing 
significantly  improved  sensitivity  critical  to  obtaining  low  energy  yields.  The  ion  source  and  QCM  are 
housed  within  a  0.125  m*  stainless  steel  vacuum  chamber  (43  cm  ID  x  76  cm  long  main  section), 
equipped  with  a  1500  liter/s  cryogenic  pump  (CTI-8).  The  chamber  base  pressure  is  5x10  *  Pa  giving  a 
working  pressure  of  approximately  0.6  to  1.3x10'^  Pa.  The  DC  four-grid  ion  source  is  specially 
designed  for  low  energy  operation  and  is  described  below  in  more  detail.  A  rotatable  target-mount  is 
positioned  23  cm  downstream  of  the  ion  source.  Both  QCM  and  target  rotation  are  performed  using 
25000  step  stepper  motors  controlled  using  a  motion  control  system  (Parker  CompuMotor).  A  personal 
computer  (LabView)  is  used  for  controlling  the  QCM  and  target  rotation  and  for  data  logging.  Detailed 
discussion  of  the  QCM  sensor  is  provided  in  Section  IID. 
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Figure  1.  Schematic  diagram  of  the  experimental  set-up. 


B.  Four-Grid  Ion  Source 

For  measurements  at  low  ion  energies  a  DC  ion  source  with  a  dual  thoriated  tungsten  fdament  as  a 
discharge  cathode  and  thoriated  tungsten  filament  as  a  neutralizer  is  used’^.  The  four -grid  ion  optics 
system,  designed  using  our  in-house  ffx  code"^'*  and  fabricated  in-house,  enables  collimated  beams  at 
the  low  ion  energies  of  interest  to  this  work  (30-350  eV).  The  grids  have  diameter  of  3  cm  and  are 
manufactured  from  Poco  graphite.  For  higher  energy  (250-1500  eV)  measurements  an  alternative  two- 
grid  ion  optics  system  with  2  cm  aperture,  also  manufactured  from  Poco  graphite,  can  be  used.  In  the 
energy  range  of  250-350  eV,  both  ion  optics  systems  can  be  used,  and  the  measurement  results  are 
consistent  between  two  systems.  The  ion  source  can  be  operated  on  different  noble  gases. 

The  four-grid  ion  optics  includes  a  screen  grid  (the  innermost),  an  acceleration  grid,  a  focusing  grid, 
and  a  deceleration  grid  (the  outermost).  The  screen  grid  is  used  to  control  the  beam  energy  and  is 
maintained  positive  to  ground,  while  the  acceleration  grid  draws  the  ion  beam  out  of  the  source  and  is 
maintained  at  negative  550-625  V  relative  to  ground.  The  focus  grid  collimates  the  beam,  and  is  held 
negative  to  ground,  generally  at  -125  V.  Finally,  the  deceleration  grid  helps  prevent  electron 
backstreaming,  and  is  generally  maintained  slightly  negative  (~2  V)  relative  to  ground.  A  typical  xenon 
flow  rate  is  0.5  seem,  typical  discharge  voltage  is  20  to  40  volts,  and  typical  neutralizer  current  is  150% 
of  the  beam  current. 
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During  sputter  measurements  the  ion  current  leaving  the  source  is  recorded.  Figure  2  shows  an 
example  of  extracted  total  beam  current  as  a  function  of  the  beam  voltage.  Monitoring  the  ion  source 
current  confirms  source  operation,  but  sputter  yields  are  determined  from  measurements  of  the  current 
reaching  the  target  (performed  before  and  after  the  sputter  tests).  When  a  conductive  target  is  being 
sputtered,  it  is  possible  to  directly  measure  the  current  reaching  the  target.  This  is  particularly  important 
for  high  incidence  angles  because,  due  to  beam  divergence,  some  fraction  of  the  ion  beam  does  not 
reach  the  target.  To  measure  the  target  current,  the  target  is  biased  negative  30  V  relative  to  the 
reference  vacuum  chamber  ground  to  repel  electrons.  For  non-conductive  targets,  the  target  current  can 
be  measured  by  replacing  the  target  with  a  metal  target  of  the  same  size.  Determination  of  the  sputter 
yield  requires  the  current  of  energetic  particles  (ions  and  fast  neutrals)  incident  on  the  sputtering  target. 
As  in  our  past  work’®  '^,  we  make  corrections  for  charge-exchange  and  scattering.  The  charge-exchange 
beam  generates  fast  neutrals  which,  depending  on  scattering  angles,  may  bombard  the  target.  The 
resulting  correction  is  to  multiply  the  measured  target  current  (due  to  ions)  by  1 .1  ±  0.05. 


Figure  2.  Beam  eurrent  vs.  beam  voltage  for  four-grid  iou  souree  opties. 


Ion  source  characterization  was  performed  in  order  to  fully  determine  the  bombardment 
conditions  (ion  energy  and  species  composition)  with  high  degree  of  accuracy.  A  collimated  Faraday 
probe  was  used  to  measure  the  beam  profile,  a  Retarding  Potential  Analyzer  (RPA)  was  used  to 
measure  the  energy  distribution  of  the  ions,  and  an  ExB  probe  was  used  to  measure  the  ratio  of  doubly 
charge  ions  to  singly  charge  ions. 


Ion  beam  profile 

Ion  beam  profiles  were  measured  using  a  guarded  Faraday  probe  attached  to  a  rotating  QCM 
mount.  The  beam  full  width  at  half  maximum  (FWHM)  angles  found  from  the  beam  profiles,  are 
shown  in  Fig.  3.  for  several  beam  voltages.  The  beam  divergence  increases  as  the  beam  energy  is 
reduced. 


Beam  voltage 

Figure  3.  Beam  eurreut  vs.  beam  voltage  for  four-grid  iou  souree  opties. 

Ion  energy  distribution 

A  large  area  four-grid  retarding  potential  analyzer  (RPA)  with  an  aperture  diameter  of  10.28 
cm  was  used  to  measure  the  ion  energy  distribution"^^.  The  large  diameter  increases  the  sensitivity  of 
the  probe,  enabling  measurements  in  low  density  plasma.  For  the  measurements  described  herein,  a 
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stainless  steel  mask  with  a  12.7  mm  diameter  opening  was  placed  on  the  RPA.  The  body  of  the  RPA 
was  made  of  aluminum  and  is  coated  with  Alodyne  to  prevent  surface  arcing.  Several  ion  energy 
distribution  functions  (lEDFs)  measured  using  the  RPA  are  shown  in  Fig.  4.  The  peak  of  the  lEDF  is 
located  within  2  V  of  the  beam  voltage  for  most  operating  conditions,  and  the  FWITM  is  6-19  volts, 
increasing  with  beam  energy. 


E,  eV 

Figure  4.  Four  differeut  ion  energy  distributions  for  beam  voltages  of  30, 80,  150  and  250 

V. 


Ion  species  measurements 

The  RPA  measurements  presented  above  do  not  discriminate  particles  by  the  charge  state.  To 
determine  the  ion  species  composition,  specifically,  the  ratio  of  doubly  charged  ions  to  singly  charged 
ions,  an  ExB  probe  was  used.  An  ExB  probe,  also  known  as  Wien  fdter,  is  a  mass  spectroscopic  device 
that  separates  charged  particles  according  to  their  velocity.  The  principle  of  the  probe  is  based  on  the 
use  of  uniform  crossed  electric  and  magnetic  fields,  which  are  perpendicular  to  each  other  and  to  the 
particle  velocity  vector.  A  schematic  of  the  probe  is  shown  in  Fig.  5.  Ions  passing  through  a  collimator 
(5  cm  long  stainless  steel  tube  with  0.5  mm  holes  on  both  sides)  (1)  enter  a  separation  region  (2),  where 
the  electric  (vertical)  and  magnetic  (perpendicular  to  the  picture  plane)  fields  act  on  it  with  Lorentz 
force.  When  a  certain  combination  of  electric  and  magnetic  fields  is  applied,  ions  with  a  corresponding 
velocity  will  pass  through  the  separation  space  and  reach  the  collector  (4).  The  suppressor  (3)  is  used  to 
suppress  secondary  electron  emission  from  the  collector  and  serves  to  improve  the  accuracy  of  the 
collector  current  measurements.  During  the  measurements,  the  voltage  (the  electric  field)  applied 
between  the  plates  in  the  separation  region  is  swept,  while  the  collector  current  is  recorded.  The 
dependence  of  collector  current  on  voltage  allows  determination  of  the  relative  number  densities  of 
different  ion  species,  namely,  singly  and  doubly  charged  xenon  ions.  In  this  work,  no  triple  charged 
ions  were  detected.  An  example  of  ExB  probe  trace  is  shown  in  Fig.  6. 


For  the  four-grid  ion  source  at  beam  energies  between  80  eV  and  250  eV  the  number  density  ratio  of 
doubly  charged  to  singly  charged  ions  is  at  most  5%,  while  for  beam  energies  below  80  eV  the  number 
density  of  doubles  is  below  the  detection  limit  (1%).  For  a  given  beam  energy,  the  doubly  charged  ions 
possess  twice  the  energy  of  their  singly  charged  counterparts.  As  shown  in  section  III,  the  yield  versus 
energy  dependence  for  the  range  of  energies  studied  is  close  to  linear,  so  that  doubly  charged  ions  have 
roughly  twice  higher  sputter  yields.  In  the  beam  current  measurement,  each  double  charged  ion  is 
counted  as  two  single  charged  ions;  therefore,  owing  to  the  relatively  small  fraction  of  double  charged 
ions,  their  effect  on  the  measured  sputter  yields  is  negligible.  To  prevent  the  formation  of  double 
charged  ions  during  low  ion  energy  measurements  (where  the  effect  on  measured  yields  could  be  more 
significant),  the  discharge  voltage  is  reduced  below  20  V  (i.e.  to  a  value  less  than  the  second  ionization 
potential  for  Xe). 
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Figure  6.  ExB  probe  traee  for  an  ion  beam  energy  of  250  eV. 


C.  Surfaee  Charging  Effeets 

Our  past  measurements  of  insulator  materials  such  as  boron  nitride  (BN)  have  shown  the 
effects  of  surface  charging  and  the  importance  of  appropriate  neutralization’ Similar  surface 
charging  effects  have  been  observed  by  Zhang  et  al."’®  and  Nikiporetz  et  al/^.  In  order  to  neutralize  the 
surface  charge,  a  plasma  bridge  neutralizer  (PBN)  has  been  placed  in  the  chamber  close  to  the  target. 
As  an  example,  Fig.  7  shows  the  dependence  of  the  (apparent)  sputter  yield  of  FIBC  grade  boron  nitride 
on  beam  current,  before  and  after  PBN  installation. 


Figure  7.  Sputter  yield  versus  beam  current  with  and  without  PBN  installation. 

As  can  be  seen  in  Fig.  7,  in  the  absence  of  surface  neutralization  there  is  a  reduction  in  sputter  yield 
with  increasing  beam  current  owing  to  charge  build-up  on  the  BN  surface.  For  higher  beam  currents, 
more  positive  charge  is  accumulated  on  the  target  surface  which  repels  beam  ions,  thus  reducing  the 
measured  sputter  yield.  Measurements  of  the  sputter  yield  versus  beam  current  with  the  PBN  show  a 
significantly  reduced  dependence  on  the  beam  current  since  the  PBN  current  neutralizes  any  surface 
charge  on  the  target  surface.  The  measurements  in  Fig.  7  were  performed  for  500  eV  xenon  ions  at 
normal  incidence  with  the  QCM  at  40“  (polar  angle)  from  the  surface  normal.  Operating  conditions  of 
the  PBN  included  an  emission  current  of  10-20  mA  and  a  Xe  mass  flow  rate  of  0.5  seem.  The  PBN 
was  biased  negatively  relative  to  ground  potential. 

D.  QCM  Sensor  and  Temperatnre  Control 

In  deposition  mode,  the  QCM  allows  determination  of  differential  sputter  yields  through 
measurement  of  mass  accumulation  (of  sputtered  particles)  on  its  surface.  For  condensable 
components,  sticking  coefficients  are  assumed  to  be  unity.  Sticking  coefficients  for  “new  layers”  and 
very  thin  layers  (on  the  order  of  several  Angstrom)  may  be  less  than  unity,  but  once  a  sufficient  layer 
thickness  of  a  given  material  has  accumulated,  sticking  coefficients  for  condensables  are  generally 
unity'’®.  For  the  multi-component  materials,  such  as  BN,  the  sputtered  particles  may  consist  of  a  mix  of 
condensable  and  non-condensable  components  and  data  must  be  interpreted  accordingly. 
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We  use  a  deposition  controller  (Sigma  Instruments  SQC-339)  that  reads  the  crystal  frequency  to 
0.001  Hz  and  an  RC-cut  quartz  crystal  as  opposed  to  the  more  conventional  AC -cut  crystal.  The  RC-cut 
crystal  (Tangidyne  Corporation)  is  very  accurate  for  deposition  of  thin  films  owing  to  increased 
sensitivity  that  is  achieved  hy  adjusting  the  stress  coefficients  of  the  quartz  plate  using  advanced 
fabrication  methods.  A  comparison  of  measured  QCM  mass  deposition  at  fixed  sputtering  conditions 
using  the  AC -cut  and  RC-cut  crystals  is  presented  in  Fig.  8.  The  ordinate  axis  is  deposited  mass,  which 
(for  fixed  sputtering  conditions)  should  increase  linearly  against  the  charge  of  ions  sputtering  the  target. 
The  data  set  measured  using  the  AC-cut  QCM  has  a  relatively  high  noise  level,  which  significantly 
reduces  measurement  sensitivity,  especially  at  low  ion  energies.  For  most  materials  of  interest,  the  RC- 
cut  crystal  enables  measurement  of  differential  sputter  yield  at  ion  energies  as  low  as  30-60  eV. 


Figure  8.  Mass  change  with  AC-cut  (left)  and  RC-cut  (right)  QCM  crystals.  Measurement 
conditions:  xenon  ions  on  Quartz  at  350  eV,  normal  incidence,  QCM  positioned  at  polar  angle  of 

20°. 


The  quartz  crystal  resonance  frequency  is  extremely  sensitive  to  temperature  so  that  for  accurate 
measurements  the  QCM  must  be  maintained  at  constant  temperature.  For  this  reason  cooling  water  is 
circulated  through  the  stainless  steel  body  of  the  QCM  housing.  Fig.  9  shows  the  importance  of 
temperature  stability,  including  temperature  induced  rippling  in  apparent  QCM  mass  change  at  fixed 
sputtering  conditions.  Again,  the  ordinate  axis  is  the  mass  accumulation,  the  slope  of  which  gives  the 
sputter  yield.  With  lesser  temperature  stability  (±0.1  K),  short  term  rippling  is  observed  and  is 
correlated  with  temperature  variation  in  the  QCM  cooling  loop.  The  effect  of  the  rippling  is  to  degrade 
our  ability  to  determine  the  slope,  thereby  degrading  overall  measurement  accuracy  and  sensitivity.  To 
remedy  this  behavior  we  have  upgraded  to  a  digital  temperature  controller  (Polyscience  9002)  with 
reffigerating/heating  circulator  that  allows  temperature  control  to  ±0.01  K.  The  right  of  Fig.  9  shows  a 
similar  measurement  with  the  improved  temperature  controller  in  which  the  rippling  is  significantly 
reduced  leading  to  improved  sensitivity  by  a  factor  of  approximately  3. 


Charge,  C 


Figure  9.  Mass  change  for  temperature  controllers  with  +0.1  K  stability  (left)  and  ±0.01  K 
stability  (right).  Measurement  conditions:  xenon  ions  on  Kapton  at  250  eV,  incidence  angle  of 
45°,  QCM  positioned  at  polar  angle  of  20°. 
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As  the  crystal  is  moved  to  different  positions  during  the  measurement,  its  (incoming)  heat  flux  varies 
due  to  changes  in  the  relative  position  of  the  QCM  and  the  heat  sources  in  the  system  (ion  source  and 
PBN).  Therefore,  although  the  temperature  of  the  water  stays  constant,  the  actual  crystal  temperature  is 
different  at  different  locations  (polar  angles).  A  K-type  thermocouple  wire  embedded  in  a  copper 
holder  silver-soldered  to  the  back  of  the  QCM  crystal  holder  is  used  to  monitor  QCM  temperature. 
When  the  QCM  is  moved  to  each  measurement  position,  the  QCM  temperature  is  monitored,  and 
sputter  yield  measurements  are  started  only  after  the  temperature  of  the  crystal  has  stabilized. 


III.  Measurement  Proeedures  and  Data  Analysis 
A.  Definition  of  Angles 

The  angles  used  to  describe  the  direction  of  ion  incidence  and  the  ejections  angles  of  sputtered 
particles  are  shown  in  Fig.  1.  We  define  as  follows;  P  is  the  incidence  angle  of  bombarding  ions 
measured  relative  to  the  surface  normal  (/5=0  for  normal  incidence),  a  is  the  ejection  polar  angle  of 
sputtered  atoms  measured  relative  to  the  surface  normal,  and  ^is  the  ejection  azimuthal  angle  of  the 
sputtered  atoms  measured  in  the  plane  of  the  target  surface  (defined  so  that  ^0  is  in  the  forward 
sputter  direction  i.e.  in  the  forward  direction  of  the  plane  containing  the  surface  normal  and  the  incident 
ion  directions). 


B.  QCM  Measurements  and  Data  Analysis 

Sauerbrey's  equation  is  used  to  calculate  the  mass  accumulation  rate  from  the  QCM  frequency 
change"^^.  A  total  of  34-36  positions  above  the  target  are  typically  sampled.  At  a  given  measurement 
point  the  volumetric  differential  sputter  yield,  y{a,  in  units  of  mmVC/sr,  is  determined  using; 


Ri^cc,  <jR)TQQM 

PJ  B,avg^s 


(1) 


where  R(a,  (j>)  is  the  measured  mass  accumulation  rate  (found  from  a  slope  of  the  mass  accumulation 
rate,  see  Figs.  8-9),  p  is  the  density  of  target  material,  JB,avg  (C/s)  is  the  time -averaged  current  of 
bombarding  particles  (ions  and  energetic  neutrals)  incident  on  the  target,  r^jcm  is  the  distance  from  the 
target  center  to  the  QCM  (17.4  cm),  and^.,  is  the  QCM  sensor  area  (0.535  cm^).  The  Cfamtity A Jr^^m 
corresponds  to  the  solid  angle  that  the  QCM  sensor  subtends  while  R(a,(f>)/pJB.avg  corresponds  to  the 
volume  of  sputtered  material  per  bombarding  charge.  It  is  important  to  emphasize  that  the  directly 
measured  quantity  is  the  mass  buildup  of  condensable  particles  on  the  QCM  and  the  volumetric 
differential  sputter  yield  should  be  considered  in  this  way.  (In  fact  the  “volume”  may  not  really 
correspond  to  any  physically  observed  volume  since  it  corresponds  to  the  equivalent  volume  due  to  the 
mass  of  the  deposited  condensable  material  if  one  uses  the  density  of  the  full  target  material,  of  course 
conversion  to  the  deposited  mass  yield  simply  requires  removing  the  density  from  Eq.  (1)). 

In  principle,  because  of  the  finite  size  of  the  QCM  crystal  and  beam  spot  on  the  target,  each  of  our 
measurements  (i.e.  QCM  positions)  corresponds  to  a  (small)  range  of  polar  and  azimuthal  angles 
joining  the  target  and  QCM.  A  simple  simulation  has  been  performed  to  show  that  for  our  chamber 
geometry  and  sputter  conditions  such  effects  are  negligible  (worst  case  of  5%  error),  so  that  the  target 
and  QCM  can  be  treated  as  points  (not  planes)  joined  by  a  single  vector®.  We  also  note  that  the 
viewing  angle  of  the  QCM  is  approximately  5.4  steradians  (cone  apex  angle  of  165°)  which 
corresponds  to  viewing  area  on  the  target  that  is  larger  than  the  irradiated  area  (also  a  necessary 
requirement  for  Eq.  (1)).  Experimental  uncertainty  and  error  bars  on  measured  differential  sputter 
yields  are  found  by  estimating  individual  contributions  and  combining  them  to  acquire  a  typical  value 
of±  16%. 

As  mentioned  above,  the  accuracy  and  sensitivity  of  QCM  measurements  is  limited  by  temperature 
drift  and  fluctuations,  which  cause  drift  and  fluctuation  of  the  QCM  signals.  The  typical  drift  rate  of  the 
QCM  signal  was  measured  when  the  QCM  temperature  controller  and  the  ion  beam  were  running  but 
no  sputtered  material  was  deposited  on  the  QCM.  The  fluctuations  of  the  QCM  signal  due  to  the 
thermal  perturbations  are  on  the  order  of  IxlO  "  g/s.  From  Eq.(l)  it  is  obvious  that,  depending  on  the 
density  of  the  material  and  the  beam  current,  this  fluctuation  can  correspond  to  different  error  margin 
and  sensitivity  threshold  in  terms  of  the  sputter  yield.  Using  gold  as  an  example  of  high  density 
material  (19.3  g/cm^)  and  with  a  high  current  value  (3  mA),  we  obtain  an  error  margin  of  IxlO'*'  g/C/sr, 
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which  is  about  two  orders  of  magnitude  lower  then  the  actual  signal  due  to  the  sputtering.  In  this  case 
thermal  drift  is  insignificant  and  the  accuracy  is  determined  by  the  other  factors.  Using  BN  as  an 
example  of  low  density  material  (2  g/  cm^),  and  with  a  low  current  value  (1mA),  we  obtain  a  value  of 
3x10'^  g/C/sr,  which  in  this  case  will  determine  the  sensitivity  threshold  of  sputter  yield  measurements. 
The  rated  sensitivity  of  the  QCM  is  least  one  order  of  magnitude  better  that  the  sensitivities  given 
above  and  work  is  underway  aimed  at  achieving  better  thermal  stability  of  the  QCM;  however,  the 
present  levels  are  adequate  for  many  measurement  needs,  even  for  low  energy  sputtering. 

C.  Modified  Zhang  Equations 

Analysis  and  fitting  of  differential  sputter  yield  profiles  requires  appropriate  functional  forms.  At 
our  conditions,  stopping  is  predominantly  due  to  elastic  (nuclear)  collisions  and  is  generally  in  the 
linear  cascade  regime  (emitted  particles  are  secondary  or  higher  generation  recoils)  or  single  knock-on 
regime  (emitted  particles  are  primary  recoils)’.  A  classical  theory  for  the  linear  cascade  regime  was 
originally  developed  by  Sigmund^”.  Independent  of  ion  incidence  angle,  the  original  Sigmund  theory 
predicts  sputtering  profiles  that  are  azimuthally  symmetric  and  approximately  diffuse  in  shape, 
corresponding  to  cosine-like  profiles  of  the  formy  cc  cos(ay  (n=l  for  a  diffuse  profile).  More  recent 
experimental  and  numerical  studies  show  a  range  of  profile  shapes.  For  normally  incident  ions  on 
polycrystalline  and  amorphous  targets,  cosine -like  profiles  are  generally  observed  with  increasingly 
under-cosine  shapes  as  ion  energy  is  lowered  and  increasingly  over-cosine  shapes  for  higher  ion 
energies  ■  •  ■  •  pot-  obliquely  incident  ions  at  relatively  high  ion  energy,  observed  profiles  also 

tend  to  be  azimuthally  symmetric.  However,  for  lower  ion  energies  the  measured  profiles  tend  to  be 
asymmetric  with  increased  sputtering  in  the  forward  direction  -  ■  -  •  Similar  profiles  have  been 

modeled  on  a  theoretical  basis^^'^"’. 

As  a  means  to  describe  the  measured  differential  sputter  yield  profiles  we  use  expressions  from 
Zhang"’^,  based  on  work  from  Yamamura  to  which  we  introduce  two  fit  parameters.  We  term  the 
resulting  expressions  as  Modified  Zhang  (MZ)”; 


TmZ  ~  ■ 


5(a) 


1- 


F  K 

— cos(^) 


1  \E 


A\E 


cos(/^)r(a)  +  -Msin(/^)sin(a)cos(^) 


r(a)  = 


3  sin 


in(«)'-l  ^cos(«)'(3sin(«)'-(-l) +  sin(a)^ 


sin  (or) 


2  sin  (or) 


1  -  sin  (or) 


(2a) 

(2b) 


where  yMz  is  the  differential  sputter  yield,  7  is  the  total  sputter  yield,  E  is  the  ion  energy,  A*  is  a 
characteristic  energy  describing  the  profile  shape,  and  the  angles  are  as  defined  above.  The  approach 
decouples  the  amplitude  of  the  angular  profiles  from  their  shape,  through  the  use  of  Y  and  E* 
respectively.  More  recent  work  by  Zhang'*’’  also  discusses  the  use  of  a  varying  energy  parameter.  In 
general,  rather  than  using  the  MZ  expressions  for  a  priori  calculation,  we  treat  Y  and  is*  as  free  fit- 
parameters  which  we  determine  from  (least-squares  fitting)  experimental  data.  Determining  the  total 
yield  in  this  way  is  roughly  equivalent  to  integrating  the  differential  sputter  yield  profiles.  Note  that 
profile  shapes  are  determined  by  the  ratio  E*/E  and  for  high  ion  energy  (E*/E  «  1)  the  MZ  expression 
reduces  to  the  diffuse  yield  (y  =  Tcos(«;)/7i).  A  least-squares  algorithm  (realized  in  MATLAB)  is  used 
to  find  the  values  of  Y  and  E*.  The  uncertainty  in  E  is  calculated  from  statistical  variation  between 
different  measurements. 


D.  Weight  Loss  Analysis 

In  some  cases  we  perform  weight-loss  measurements  as  a  complimentary  means  to  determine 
total  sputter  yields.  For  multi-component  materials  with  sputter  products  that  may  contain  a  mix  of 
condensable  and  non-condensable  species  (e.g.  BN),  comparison  of  total  yields  from  QCM  (which 
detects  only  condensables)  versus  weight-loss  (which  detects  all  species)  provides  additional  insight 
into  the  identity  of  sputter  products  (i.e.  the  fraction  that  is  condensable).  A  microgram  scale  (Mettler 
AE163)  is  used  to  weigh  the  targets  before  and  after  each  sputtering  session.  With  this  method  the 
volumetric  sputter  yield,  Y  (in  units  of  mmVC),  of  each  target  is  found  using; 


7  = 


Am 

P  B,Avg  ^ 


(3) 
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where  Am  represents  the  mass  loss  due  to  ion  beam  exposure,  t  is  the  exposure  time,  and  other  symbols 
are  defined  as  above. 

Measurement  uncertainties  for  total  sputter  yields  are  calculated  by  estimating  individual 
contributions  and  combining  them.  Uncertainties  are  as  follows;  +  15%  on  beam  current,  +  200  |ig  on 
mass  measurements,  and  ±1%  on  test  time  lengths.  The  error  on  mass  tends  to  dominate  at  low  yields 
and  our  mass  uncertainty  is  an  order  of  magnitude  larger  than  that  provided  by  the  scales  themselves 
owing  to  the  observed  reproducibility  in  our  experiment  (see  below).  The  resulting  uncertainty  in  total 
yields  varies  depending  on  measurement  conditions,  but  typical  values  are  approximately  25%.  Note 
that  this  value  is  a  random  error  resulting  from  finite  accuracy  of  the  mass,  current,  and  time 
measurements. 

One  factor  that  can  affect  the  weight  loss  measurement  accuracy  is  moisture  absorption  by  the 
target  materials,  for  example  some  grades  of  boron  nitride.  When  exposed  to  atmosphere,  the  targets 
start  absorbing  moisture,  which  increases  their  weight.  In  order  to  minimize  this  effect,  the  vacuum 
chamber  is  vented  with  pure  nitrogen  and  the  time  of  atmospheric  exposure  before  the  weight 
measurement  is  minimized.  Another  factor  that  might  affect  accuracy  of  the  weight  loss  measurements 
is  xenon  ion  implantation  into  the  target.  If  a  significant  number  of  xenon  ions  are  implanted  in  the 
target,  their  effect  would  be  to  reduce  the  measured  mass  loss  and  the  corresponding  sputter  yield.  In 
order  to  check  the  influence  of  ion  implantation  in  one  experiment.  X-ray  photoelectron  spectrometry 
(XPS)  analysis  was  performed  on  a  BN  target  (PHI  5800  X-ray  photoelectron  spectrometer).  XPS 
measurement  confirmed  that  some  xenon  implantation  did  occur  in  a  thin  near-surface  layer  with 
thickness  not  exceeding  a  few  nanometers.  Simple  calculation  shows  the  mass  contribution  from 
implanted  xenon  to  be  negligible  compared  to  the  mass  change  due  to  sputtered  material. 

E.  Measurement  Proeedure 

Before  the  measurements,  targets  are  deliberately  pre -sputtered  to  remove  any  surface  layer  that 
could  have  been  formed  by  contamination  or  mechanical  modification  during  the  manufacturing 
process.  Pre -sputtered  surfaces  also  better  represent  the  conditions  found  in  long-duration  electric 
propulsion  operating  applications.  An  order -of-magnitude  estimate  for  the  typical  dose  of  incident  ions 
on  a  target  prior  to  testing  is  approximately  10*^  ions/cm^  (corresponding  to  thickness  of  several 
microns).  Target  contamination  effects  in  the  vacuum  chamber  are  estimated  to  be  negligible,  since  for 
typical  conditions  the  flux  of  ions  incident  on  the  target  is  approximately  10  times  higher  than  the  flux 
of  nitrogen  (the  maj  or  contaminant)  to  the  target*. 

Before  placement  into  the  vacuum  chamber,  targets  are  placed  on  the  scales  described  above.  Each 
target  is  weighed  5  times  to  ensure  stability,  and  the  average  value  of  the  measured  mass  is  taken. 
Targets  are  then  inserted  into  the  vacuum  chamber  and  the  sputter  measurement  is  conducted.  Test 
durations  are  fixed  such  that  the  mass  change  of  the  test  is  large  compared  to  the  precision  of  mass 
measurements,  and  such  that  the  QCM  has  time  (at  each  position)  to  sufficiently  stabilize  relative  to 
thermal  and  background  noise.  Test  times  vary  from  less  than  an  hour  for  higher  sputter  yield  materials 
(e.g.  Mo,  Au)  and  higher  energy  tests  (e.g.  350  eV)  to  as  long  as  15  hours  for  lower  sputter  yield 
materials  (e.g.  BN)  and  lower  energy  tests  (e.g.  50  eV).  After  each  test,  the  targets  are  weighed  again 
and  the  difference  in  masses  is  used  for  the  weight-loss  measurement. 


IV.  Validation  measurements 

In  this  section  we  present  validation  measurements  demonstrating  system  sensitivity  and  accuracy. 
The  sputter  yield  of  Xe  ions  on  Au  was  measured  for  the  ion  energy  range  30-350  eV  and  at  incidence 
angles  of  0-67.5°.  The  total  yields  from  weight  loss  and  QCM  (found  from  the  Y  parameter  of  the  best- 
fit  MZ  profile)  are  self-consistent  to  within  experimental  uncertainty.  Fig.  10  presents  the  total  yield 
versus  ion  energy  as  well  as  literature  data  for  comparison**’**.  The  energy  dependence  is  generally  in 
good  agreement  with  the  Yamamura  and  Tawara  empiric  model**.  It  should  be  noted  that  the  threshold 
sputter  energy  for  Xe  on  Au  predicted  by  Yamamura  and  Tawara  is  19  eV,  which  is  also  in  broad 
agreement  with  our  measurements  -  the  lowest  energy  where  sputtering  was  detected  is  30  eV.  Both 
our  data  and  the  Yamamura  and  Tawara  predictions  are  consistently  higher  than  the  Wehner  and 
Rosenberg  data**,  though  similar  trends  are  observed  in  all  cases. 
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Figure  10.  Total  sputter  yield  of  gold  due  to  xeuon  Ions  at  normal  ineidenee. 

The  dependence  of  the  total  sputter  yield  on  ion  incidence  angle  for  xenon  ions  on  gold  is  given  in 
Fig.  11.  The  sputter  yield  increases  up  to  incidence  angle  of  45°,  and  then  drops,  which  is  roughly 
characteristic  for  most  materials.  A  prediction  from  the  Yamamura  and  Tawara  model^^  is  also  given 
and  is  in  reasonable  agreement  with  our  measurements,  though  the  model  predicts  the  peak  sputter 
yield  to  be  at  a  slightly  higher  incidence  angle. 


Figure  11.  Total  sputter  yield  of  gold  due  to  xenon  Ions  versus  ineidenee  angle. 

Examples  of  differential  sputter  yield  profiles  due  to  xenon  ions  on  gold  are  shown  in  Figs.  12  and 
13.  Figure  12  shows  differential  sputter  yield  profiles  for  350  eV  ions  at  normal  incidence,  and  for  250 
eV  ions  at  15°  incidence  ((|)=0/180°  plane).  Experimental  data  and  the  corresponding  modified  Zhang 
fits  are  shown,  as  well  as  diffuse  profiles  (with  the  same  total  sputter  yields)  for  comparison.  For  both 
conditions  shown  in  Fig.  12,  as  well  as  for  all  other  conditions  studied  in  this  work,  the  differential 
sputter  yield  profiles  were  significantly  non-difflise,  emphasizing  the  importance  of  measuring  the  full 
angular  profiles  as  opposed  to  assuming  a  diffuse  shape.  For  normally  incident  ions,  the  observed 
profiles  for  the  low  (<500  eV)  energy  ions  are  typically  “under-cosine”,  meaning  reduced  sputtering  in 
the  normal  direction  as  compared  to  a  diffuse  profile,  as  in  Fig.  12a.  For  non-normal  incidence,  the 
profiles  are  typically  azimuthally  asymmetric  with  increased  sputtering  in  the  forward  ((|)=0)  direction. 
Fig.  13  shows  the  hemispheric  representation  of  the  differential  sputter  yield  from  the  best-fit  modified 
Zhang  equation  for  the  same  two  cases. 
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a,  degrees 

Figure  12.  Differential  sputter  yield  profiles  xenon  ions  on  gold,  (a)  350  eV  ions  at  normal 
ineidenee.  (b)  250  eV  ions  at  15“  Ineldenee,  (|f0/180“  plane. 


Figure  13.  Zhang-fit  hemisphere  for  xenon  ions  on  gold,  (left)  350  eV  ions  at  normal  ineidenee. 
(right)  250  eV  ions  at  15“  ineidenee,  (|)=0/180“  plane.  Ion  ineidenee  dlreetlon  is  shown  with  arrow. 

Differential  sputter  yields  are  in  units  of  mm^/C/sr. 

The  system  has  been  used  to  measure  total  and  differential  sputter  yields  of  a  variety 
materials,  ineluding  metals  (e.g.  molybdenum,  gold),  ceramics  (e.g.  different  grades  of  boron  nitride), 
polymers  (e.g.  Kapton),  and  other  multi-component  materials  (e.g.  quartz).  Examples  of  hemispheric 
representation  of  the  differential  sputter  yields  from  the  best-fit  modified  Zhang  equation  for  quartz  and 
Kapton  are  presented  in  Fig.  14.  For  relatively  low  ion  energies  (<~500  eV),  all  studied  materials 
display  a  non-diffuse  differential  sputter  yield  profiles.  The  energy  dependence  of  the  total  sputter 
yields  of  several  materials  is  given  in  Fig.  15,  demonstrating  the  wide  sensitivity  range  of  the 
measurement  system.  Note  that  the  yield  values  for  the  multi-component  materials  in  Fig.  15 
correspond  to  QCM-measured  yield  of  condensable  components. 


Figure  14.  Zhang-fit  hemisphere,  (left)  350  eV  xenon  ions  on  quartz  at  30°  ineidenee.  (right)  500 
eV  xenon  ions  on  Kapton  at  45“  Ineidenee.  Ion  ineidenee  direetlon  is  shown  with  arrow. 
Differential  sputter  yields  are  in  units  of  mm^/C/sr. 
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V.  Conclusions 

We  have  deseribed  a  QCM  based  system  used  for  differential  sputter  yield  measurements.  For  low 
ion  energies,  the  system  uses  a  four  grid  ion  source.  Source  characterizations  were  presented  and  show 
monoenergetic  beams  with  low  fractions  of  multiply  charged  ions.  The  QCM  employs  a  high- 
sensitivity  RC-cut  crystal  to  measure  the  amount  of  material  sputtered  in  a  given  angular  direction. 
Scanning  the  QCM  over  the  hemisphere  allows  measurement  of  the  full  differential  sputter  yield 
profile.  Temperature  stability  is  found  to  be  critical  and  a  controller  with  ±0.01  K  accuracy  is  used  to 
stabilize  QCM  temperature.  A  plasma  bridge  neutralizer  is  used  for  non-conductive  targets  to  prevent 
positive  charge  buildup  on  the  target  surface.  Validation  measurements  for  Xe  ions  on  Au  are  presented 
for  ion  energies  of  30-350eV  and  incidence  angles  of  0-67.5“.  Total  sputter  yields  from  integrating  the 
QCM  differential  sputter  yields  are  consistent  with  weight  loss  measurements  and  in  agreement  with 
published  values.  Examples  of  results  from  other  materials  are  also  presented. 

Despite  the  need  for  low  energy  (<100  eV)  differential  sputter  yields  in  a  range  of  applications, 
there  is  currently  a  lack  of  suitable  measurement  systems  for  this  purpose.  The  high-sensitivity  and 
directly  quantitative  nature  of  the  QCM  based  system  presented  here  contribute  to  filling  this  gap. 
Differential  sputter  yield  profiles  are  presented,  and  for  the  conditions  studied,  under-cosine  profiles 
are  observed  at  normal  incidence,  while  profiles  with  increased  forward  sputtering  are  observed  for 
non-normal  incidence.  In  all  cases,  non-diffiise  profiles  are  observed  showing  the  need  to  directly 
measure  the  angular  dependence  of  the  sputter  yields. 
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